Irwin et al. BMC Plant Biology 2012, 12:21 
http://www.biomedcentral.eom/1 471-2229/1 2/21 



Plant Biology 



RESEARCH ARTICLE Open Access 



Functional alleles of the flowering time regulator 
FRIGID A in the Brassica oleracea genome 

Judith A Irwin 1 , Clare Lister 1 , Eleni Soumpourou 1 , Yanwen Zhang 1 , Elaine C Howell 2 , Graham Teakle 3 and 
Caroline Dean 1 " 



Abstract 

Background: Plants adopt different reproductive strategies as an adaptation to growth in a range of climates. In 
Arabidopsis tholiona FRIGIDA {FRI) confers a vernalization requirement and thus winter annual habit by increasing 
the expression of the MADS box transcriptional repressor FLOWERING LOCUS C {FLQ. Variation at FRI plays a major 
role in A. tholiono life history strategy, as independent loss-of-function alleles that result in a rapid-cycling habit in 
different accessions, appear to have evolved many times. The aim of this study was to identify and characterize 
orthologues of FRI in Brassica oleracea. 

Results: We describe the characterization of FRI from Brassica oleracea and identify the two B. oleracea FRI 
orthologues (BolCFRI.a and BolCFRI.b). These show extensive amino acid conservation in the central and C-terminal 
regions to FRI from other Brassicaceae, including A. thaliana, but have a diverged N-terminus. The genes map to 
two of the three regions of B. oleracea chromosomes syntenic to part of A thaliana chromosome 5 suggesting 
that one of the FRI copies has been lost since the ancient triplication event that formed the B. oleracea genome. 
This genomic position is not syntenic with FRI in A thaliana and comparative analysis revealed a recombination 
event within the A thaliana FRI promoter. This relocated A thaliana FRI to chromosome 4, very close to the 
nucleolar organizer region, leaving a fragment of FRI in the syntenic location on A thaliana chromosome 5. Our 
data show this rearrangement occurred after the divergence from A lyrata. We explored the allelic variation at 
BolCFRI.a within cultivated B. oleracea germplasm and identified two major alleles, which appear equally functional 
both to each other and A. thaliana FRI, when expressed as fusions in A. thaliana. 

Conclusions: We identify the two Brassica oleracea FRI genes, one of which we show through A. thaliana 
complementation experiments is functional, and show their genomic location is not syntenic with A. thaliana FRI 
due to an ancient recombination event. This has complicated previous association analyses of FRI with variation in 
life history strategy in the Brassica genus. 
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Background 

The switch to reproductive development is a fundamen- 
tal process in the plant life cycle. The molecular mechan- 
isms underlying this developmental transition have been 
extensively studied in Arabidopsis thaliana. An inte- 
grated network of environmentally responsive genetic 
pathways converge on a common set of targets to quanti- 
tatively regulate the genes required to switch the apical 
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meristem from a vegetative to a floral state [1-3]. One 
important environmental cue is prolonged cold, which 
accelerates flowering in a process termed vernalization 
and aligns pollination and seed set with the favourable 
conditions of spring. Variation in requirement for verna- 
lization exists in many plant species and this influences 
life history strategy with plants requiring vernalization 
adopting a perennial, biennial or winter annual habit in 
contrast to summer annuals, which flower in the first 
growing season. This is in contrast to other species that 
are more reliant on photoperiodic signals or endogenous 
cues e.g. rice [4]. The significant fitness consequences of 
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flowering time variation, demonstrated in annual [5,6] 
and perennial plants [7], have most likely contributed to 
the evolution of the extensive variability in flowering 
time control. Flowering also influences the pattern of 
growth throughout the seasons and affects many agro- 
nomic characters including the quantity and quality of 
crop production. This is particularly apparent in culti- 
vated brassicas, where variation in the flowering process 
has been selected to produce a diverse array of econom- 
ically important morphological forms. 

A major determinant in the variation of vernalization 
requirement in A. thaliana is allelic variation at FRIGIDA 
{FRI) [8-11]. FRI represses flowering by promoting the 
expression of the floral repressor FLOWERING LOCUS C 
(FLC) [12,13]. Vernalization acts antagonistically to FRI 
and accelerates flowering by down-regulating FLC. A 
number of rapid-cycling variants of A. thaliana that do 
not need vernalization were found to have arisen through 
loss of function of FRI, an evolutionary step that has 
occurred multiple times [8,9,11,14]. Parallel evolution 
through allelic variation at a common target has been 
found in other organisms [15]. It was therefore interesting 
to ask whether a similar evolutionary step has occurred in 
other plant species. Many other species do show variation 
in vernalization requirement and it is an important agro- 
nomic trait in many major crops. For example, in B. olera- 
cea (horticultural brassicas) vernalization-requiring 
biennials are represented by cabbage and Brussels sprouts, 
with summer annual crops including some calabrese and 
cauliflower cultivars. Orthologues of FRI have been identi- 
fied in A. lyrata [16], Capsella species [17] and the halo- 
phyte Thellungiella halophila [18] within the Brassicaceae, 
and more broadly in Medicago truncatula, Lotus japoni- 
cus, Vitis vinifera [19], Populus balsamifera [20] and 
Oryza sativa [21]. To date natural variation in vernaliza- 
tion requirement has been associated with FRI polymorph- 
ism in A. lyrata [8] and allelic variation in one orthologue 
in Brassica napus (BnaA.FRI.a) has been associated with 
flowering time variation [22] . 

We are interested in understanding the molecular basis 
of variation in flowering time and vernalization require- 
ment in horticultural brassicas. Genetic information from 
A. thaliana can generally be applied to Brassica species 
because of their evolutionary relatedness. The Arabidop- 
sis and Brassica genera are in the same family (Brassica- 
ceae) with B. oleracea thought to have arisen from a 
triplication of an ancestral genome similar to that of Ara- 
bidopsis [23-26]. Genetic information on the control of 
flowering in Arabidopsis can be applied to Brassica spe- 
cies because of the colinearity of the Arabidopsis and 
Brassica genomes [23,27,28]. This has been used to infer 
candidate genes that might account for QTL underlying 
flowering time and other variation [22,29,30]; however, in 
some instances it can be misleading [31]. Here, we 



identify the two FRI genes in the B. oleracea genome and 
map their genomic locations. We also explored allelic 
variation at one of the FRI loci in cultivated B. oleracea 
germplasm. These new data will provide the necessary 
information to elucidate how general a role FRI plays in 
life history variation in the Brassicaceae. 

Results and discussion 

Two FRI genes are present in the Brassica oleracea 
genome 

The BoFRI genes were isolated from the JBo BAC library 
of the B. oleracea Chinese kale genotype A12DHd [32] 
through hybridization with an A. thaliana FRI genomic 
clone. From seven positive BAC clones two that showed 
distinct FRI hybridization patterns (JBo72I23 and 
JB088GI6, Figure la, b) were selected for sub-cloning. 
Analysis of these confirmed they carried different Bras- 
sica paralogues designated BolCFRI.a and BolCFRI.b 
[33] and referred to hereafter as BoFRIa and BoFRIb 
[GenBank JN191450 and JN191449]. As in other species, 
BoFRIa and BoFRIb contain three exons encoding pre- 
dicted open reading frames (ORFs) of 594 and 585 resi- 
dues respectively (Figure lc). BoFRIa contains two 
coiled-coil domains, typically involved in protein oligo- 
merisation (as predicted by COILS http://www.ch.emb- 
net.org/software/COILS_form.html [34]), very similar to 
the predicted structure of the A. thaliana FRI (AtFRI) 
[8,35]. In contrast, BoFRIb is predicted to contain only 
one coiled-coil domain in the C-terminal region as was 
found to be the case for two of the four FRI identified in 
B. napus [22]. 

AtFRI is the original member of a family of seven pro- 
teins in A. thaliana which, apart from the two predicted 
coiled-coil domains, show no homology with any other 
proteins and whose function has yet to be determined. 
Recent analysis of the FRI protein family [19] identified a 
conserved core central domain. Outside of this domain 
significant variation is observed that allows the FRI family 
to be subdivided into five distinct groups. AtFRI and its 
orthologues in other species are defined by a conserved 
region of 37 amino acids in the N-terminal region of the 
protein. The BoFRI proteins we describe here contain this 
conserved 37 amino acid region reinforcing the view they 
are FRI orthologues; however, the amino acids either side 
of this region show lower homology to AtFRI (Figure lc). 
This region includes much of the first predicted coiled- 
coil in BoFRIa. Variations in this domain in BoFRIb result 
in the loss of a predicted coiled-coil, emphasising a possi- 
ble functional significance for the amino acid polymorph- 
isms in this region. A similar degree of divergence from 
AtFRI is found in the N-terminal region of an orthologue 
of FRI isolated from the halophyte T. halophila and in 
four orthologues of FRI identified in B. napus [18,22]. By 
contrast, there is extensive amino acid conservation 
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Ji V1TAFKRR YDEL Q SHI ID YIEHflll) SHLKTHGIV -E I MVS PPPKDBSGET ATBI AC -QSPPK 

MID BFKRHYDDL QHHMD YIKNAI GS SLKFKGI IAESPSSRS — Q SPRHDaSGETMBVMTQSPPKETSETVPE - 

BED BFKRHYDDL QKHtlDD EEHflXE SKFKSHG -V -DD - -SSS — H SPEHD ASREI AT AIVC - -PPPPEEAET BPEM 
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BoFRIa IRQ S GT BE IVGflLKRS PFL VP1 « G GIVD S S TKRGHHIE BLFJ JVYTF GHEDRF SPSSLLT SFLR1 ISKE SFERflKRKflQ API! 
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AlpFRI IRMS GSHE I flGBLKRS PFLVPIH S GIVE S C TKRG1 HEE flLFJ JVYTF GMEDKF S AS SVL T SFLR1 ISKE SEERTKRKBQ S PL 
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AIIFRI RSVYA YEHL P PIISYS P GHGQRL PRQYS P S P VHGQRHPRQ YS P P THGQ QQIPFG-L QRVYRHS P SEEP YL GL SHHRS PRSH 
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Figure 1 Cloning BoFRIa and BoFRIb. (a) Southern analysis of A12 BACs identified by colony hybridisation and probed with AtFRI. Lanes 1, 4, 5, 8, 11, 

12 and 13 contain clones that show homology to FRIGIDA. BACs in lanes 2, 3, 6, 7, and 10 do not cross hybridise (b) Hindlll digest of six BACs probed 
with conserved region from exons 2 and 3 of BoFRIa. Lane 5 contains JBo72l23 from which BoFRIa was sequenced. Lane 7 contains JB088GI6. Lanes 1, 
2, 3 and 6 contain four further BACs showing the same hybridization pattern as JB088GI6. Note the intensity of the hybrization is indicative of the 
sequence divergence between BoFRIa and BoFRIb (See 1c) (c) Comparison of the protein sequences of BoFRIa and BoFRIb with other members of the 
FRI sub-family. From top to bottom they are Brassica oleracea BoFRIa, Brassica oleracea BoFRIb, Fhelliungiella halophila ThFRI, Arabidopsis thaliana AtFRI 
Arabidopsis lyrata ssp lyrata AIIFRI and Arabidopsis lyrata ssp petraea AlpFRI. The N-terminal domain containing the conserved region of 37 amino acids 
(indicated by solid bar) that defines copies of FRIGIDA from other members of the FRI superfamily [17]. The coiled-coil domains are indicated by the 
red lines. 
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between BoFRIa, BoFRIb and AtFRI in the central and C- 
terminal regions (Figure lc). Transgenic analysis of the 
functional domains of AtFRI in A. thaliana where the N 
or C terminus was deleted revealed that the N-terminal 
region was less important for function [19], perhaps 
explaining the high degree of divergence observed. 

The BoFRI genes map to regions that are non-syntenic 
with A thaliana FRI 

A genomic fragment including exon 2, intron 2 and 
exon 3 of BoFRIa (and showing a high level of conserva- 
tion in BoFRIb and AtFRI) was hybridized to mapping 
filters from two B. oleracea mapping populations: Chi- 
nese kale x calabrese (var. alboglabra x var. italica; 
A12DHdxGDDH33, [36]; Figure 2a) and cauliflower x 
Brussels sprout (var. botrytis x var. gemmifera; N x G 
[37]). RFLPs for one of the two BoFRI loci segregated in 
the A12DHdxGDDH33 mapping population that 
allowed this locus to be mapped to 39.5 cM on linkage 
group C3 of the B. oleracea genome. The locus mapping 
to C3 was identified as BoFRIa by fluorescence in situ 
hybridization (FISH) with BAC JBo72I23, from which 
BoFRIa was originally sequenced (Figure 2b). JB088GI6 
was located on the short arm of chromosome C9 by 
FISH (Figure 2c, d). Therefore, the second locus, 
BoFRIb, was on linkage group C9. Two further BACs 
showing the same restriction pattern as JB088GI6 
(Figure lb) hybridized to the same location on C9 (data 
not shown). These results confirm that the B. oleracea 
genome contains two orthologues. 

Comparative analysis of the Brassica and A. thaliana 
genomes has shown that the chromosomal regions of C3 
and C9 to which the two BoFRI loci have been mapped 
are syntenic to a region of A. thaliana chromosome 5 
and not to the top of A. thaliana chromosome 4, where 
AtFRI (At4g00650) is located [23]. This region of chro- 
mosome 5 includes a number of genes known to be 
involved in the control of flowering including FLC, FY 
and CONSTANS (CO). Several QTL studies have found 
loci for flowering time variation mapping to this genomic 
region in a number of Brassica populations including 
B. oleracea [38-41], B. rapa [42-46]; B. nigra [47,48] and 
B. napus [29,49]. The mapping we have undertaken 
reveals the proximity of BoFRI not only to BoFLC } but 
also BoFY and BoCO; other flowering time genes that 
have been mapped previously. The sequences of BoFRIa 
and BoFRIb further allow us to identify which of the four 
orthologues of AtFRI recently identified in B. napus [22] 
are the two C genome copies. The four copies of FRI 
were designated BnaA.FRI.a, and BnaX.FRIb-d. Compari- 
son of the amino acid sequences of these proteins with 
BoFRIa and BoFRIb suggest that BnaX.FRI.d is the ortho- 
logue of BoFRIa and the C genome homoeologue of 
BnaA.FRI.a. This conclusion is further supported by the 
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Figure 2 Mapping BoFRIa and BoFRIb. (a) A1 2DHdxGDDH33 
mapping population probed with conserved region from exons 2 
and 3 of BoFRIa. Two loci are identified; one monomorphic (upper 
band) and a second segregating with the two parental alleles 
(lower two bands), (b) Meiotic pachytene spread with JBo72l23 
(BoFRIa, green, white arrow) hybridizing to C3 between the 
telomere and JBo62M08 (red, yellow arrow), (c) and (d) JB088GI6 
(BoFRIb) hybridizes to the short arm of C9: (c) mitotic metaphase 
with JB088GI6, (green, arrows), BoB061G14 and 45S rDNA (red), (d) 
reprobe of (c) with JBo32J18 (green, arrows), a marker for C9. 



fact that BnaA.FRI.a was mapped to a region of A3 
homoeologous to the region of C3 where we have 
mapped BoFRIa. Comparison of the amino acid sequence 
of BnaX.FRI.c shows it to be identical to that of BoFRIb. 
BnaX.FRI.c appears most similar to BnaX.FRI.b and is 
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therefore likely to be the A genome homoeologue of 
BoFRIb in B. napus. 

A recombination event specific to the A thaliana lineage 
has relocated the FRIGIDA gene to the top of 
chromosome 4 

In A. thaliana, the AtFRI locus is located at the top of 
chromosome 4. However, it has previously been 
reported that the orthologue of FRI in A. lyrata maps to 
linkage group 8 [50,51]. This linkage group is ortholo- 
gous to the lower arm of A. thaliana chromosome 5 
[50,52,53]. Interestingly, an annotated gene model in 



this region of A. thaliana chromosome 5 (AtSgS1090) 
shows a high degree of homology to AtFRI, containing 
parts of intron 1 and exon 3 but lacking other parts of 
the coding region, thus it may be a pseudogene [50]. 
Genevestigator data suggest AtSgS1090 is expressed at 
very low levels, supporting this hypothesis [54]. Down- 
stream, in the opposite orientation, is AtSgSllOO, encod- 
ing an iron superoxide dismutase and the BoFRIa BAC 
clone contains 3' sequence showing homology to exons 
3-9 of this A. thaliana chromosome 5 gene (Figure 3a). 

Synteny has been studied extensively in the Brassica- 
ceae genomes due to its potential for gene identification 
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Figure 3 BoFRIa and BoFRIb map to regions of the B. oleracea genome that are non-syntenic with A. thaliana. (a) Comparison of BoFRIa 
genomic clone with annotated regions of A thaliana chromosomes 4 and 5 represented in the 5' to 3' orientation. Solid black rectangles 
correspond to the three exons that make up FRI. Arrows indicate the orientation of other genes, (i) BoFRIa, genomic clone, (ii) Region of 
A thaliana chromosome 5 containing an iron superoxide dismutase {FSD2, At5g5ll00) and At5W90 showing homology to FRI intron 1 and exon 
3. (iii) Region of A thaliana chromosome 4 containing AtFRI and At4g400651 showing homology to exon 9 of FSD2. (b) Proposed derivation of 
A. thaliana and B. oleracea linkage groups from ancestral karyotype (AK6 and AK8) as suggested in [53] showing the position of FRI and FLC. The 
rearrangement of ancestral blocks W and R results in orthologues of FRI and FLC being brought into close proximity in the ancestral Brassica 
genome. This region is present on chromosomes C2, C3 and C9 of B. oleracea and is located in regions of the B. oleracea genome showing 
synteny to A. thaliana chromosome 5. BoFRIa is located on C3 and BoFRIb on C9. The copy of FRI on C2 has been lost. 
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and marker development. Arabidopsis and Brassica are 
thought to have diverged about 43 Mya with a triplica- 
tion of an ancestral genome (similar to that of Arabi- 
dopsis) occurring approximately 23 Mya and giving rise 
to modern day diploid Brassica genomes [55]. A. thali- 
ana and A. lyrata are thought to have diverged around 
13 Mya, with a reduction in chromosome number, from 
the ancestral karyotype of n = 8 (as represented in A. 
lyrata) to the derived state in A. thaliana of n = 5 
[24,55,56]. The ancestral karyotype of the Brassicaceae is 
proposed to be eight chromosomes composed of 24 
conserved chromosomal blocks [57]. These blocks can 
be rearranged to model the genome structure of A. 
thaliana, A. lyrata and the modern day diploid Brassicas 
[24]. Thus the genomic composition of the nine chro- 
mosomes (C1-C9) of B. oleracea and ten chromosomes 
(A1-A10) of B. rapa can be related to both the ancestral 
karyotype and the A. thaliana genome. 

The ancestral genomic blocks QR and WX from chro- 
mosomes 6 and 8 respectively of the ancestral karyotype, 
and today represented by A. lyrata [24], have been 
recombined in the ancestral Brassica genome prior to tri- 
plication, leading to the block WR being represented 
three times in the B. rapa genome on A2, A3 and A10 
[25,26]. The paralogous regions of B. oleracea are on C2, 
C3 and C9 (Figure 3b). This rearrangement brings ortho- 
logues of FLC (block R) and FRI (block W) together on 
these chromosomes. Thus in A. lyrata and B. oleracea, 
FRI maps quite closely to VIN3 (also in block W and 
required for vernalization [58], as well as its major target 
FLC (block R; Figure 3b). BoFRIa mapping to C3 and 
BoFRIb to C9 thus represent two of the three syntenic 
regions. The third paralogue of FRI appears to have been 
lost from C2 during B. oleracea evolution; such gene loss 
is not uncommon [59]. This is in contrast to the current 
location of FRI at top of chromosome 4 in A. thaliana 
that shows homology to block O from chromosome 6 of 
the ancestral karyotype. 

The data we present here suggest that the chromoso- 
mal rearrangements that occurred during the evolution 
of the ancestral Brassicaceae genome into A. thaliana 
included a recombination/rearrangement event that 
relocated a genomic region containing AtFRI to a posi- 
tion near the distal end of the short arm of chromosome 
4, close to the nucleolar organiser region, leaving a non- 
functional remnant in the genomic position on chromo- 
some 5 that is syntenic with FRI in the other Brassica- 
ceae {A. lyrata, B. oleracea; [25,52]). 

Two common alleles exist for BoFRIa in diverse 
genotypes of B. oleracea 

The original sequences of BoFRIa and BoFRIb were 
obtained from A12DHd, one of the parents of the map- 
ping population used in [38,39]. These studies mapped 



QTL for flowering time on C3 in the region where we 
have mapped BoFRIa. We therefore sequenced 650 bp 
from exon 1 of BoFRIa from the other parent of this 
population, GDDH33 (data not shown). The GDDH33 
sequence showed two amino acid substitutions (A118V 
and Q125E) compared to A12DHd. Thus the parents of 
this population are carrying different alleles of BoFRIa 
and it is possible that variation at BoFRIa is contributing 
to variation in flowering time in this population. Single 
amino acid substitutions have also been identified in 
alleles of BnaA.FRI.a sequenced from the parents of the 
Tapidor x Ningyou7 (TN) mapping population from 
B. napus and mapped to a region underlying a QTL for 
flowering time variation [22]. 

We sequenced BoFRIa and BoFRIb for two additional 
genotypes of B. oleracea italica. The A12DHd reference 
sequence is derived from a BAC clone of a Chinese 
kale, B. oleracea alboglabra, which flowers after 8 weeks 
[38] and can be considered a rapid-cycling type, not 
requiring vernalization. We therefore selected two addi- 
tional genotypes of broccoli (Brassica oleracea italica); 
El, which has a facultative vernalization response, flow- 
ering earlier following a period of cold, but which 
matures in October/November (Autumn) of the year of 
planting and E8 which has an obligate vernalization 
requirement and matures in April/May (Spring) of the 
following year. BoFRIb is highly conserved between the 
three genotypes with only 5 amino acid substitutions 
(D6G, K20Q, Q372K and R532W in El and N105T in 
E8). The sequencing of BoFRIa in these genotypes iden- 
tified a polymorphic region in exon 1 that included two 
deletions of seven and three amino acids in E8 relative 
to El, either side of the conserved block of 37 amino 
acids that defines the FRI proteins, (Figure lc). Thirteen 
non-synonymous and 12 synonymous substitutions dif- 
ferentiate the A12DHd and El BoFRIa alleles from the 
E8 allele. We therefore designated the El and E8 
BoFRIa alleles as BoFRIa-I and BoFRIa-4 respectively 
[GenBank JN191393, JN191392]. 

We focused our subsequent analysis on BoFRIa as this 
showed most polymorphism and extended our analysis 
to include 55 genotypes from the cultivated B. oleracea 
Diversity Foundation Set developed at the University of 
Warwick (BolDFS, King et al. http://www.brassica.info/ 
resource/plants/diversity_sets.php; [60] Table 1). A 650 
bp region of BoFRIa covering the exon 1 polymorphic 
region containing the two deletions was sequenced 
(Table 1, GenBank JN191394-191448). We identified six 
BoFRIa alleles within this subset of 55 genotypes from 
the BolDFS. These can be divided into two groups; 
BoFRIa 1-3 and BoFRIa 4-6 where alleles in the second 
group include the seven amino acid and three amino 
acid deletions. The BoFRIa-1 and BoFRIa-4 alleles were 
the most common within the 55 genotypes studied. In 
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addition the two deletions identified in BoFRIa-4 always 
co-occurred and were found at high frequency together 
with a small number of non-synonymous nucleotide 
polymorphisms. The two deletions, which had been 
found in the late-flowering broccoli, were over-repre- 
sented in B. oleracea vegetable types such as Brussels 
sprouts and kohl rabi with a winter annual or biennial 
habit, usually grown for consumption of their vegetative 
rather than floral forms (Table 1, Figure 4). Interestingly 
BnaX.FRLd from the B. napus winter variety Express 
[22], which we have identified here as the C genome 
homologue of BoFRIa in B. napus has both of the dele- 
tions identified in the BoFRI 4-6 class of alleles that are 
overrepresented in brassica vegetable types with a winter 
annual or biennial habit. On closer examination BnaX. 
FRI.d was found to have the same amino acid sequence 
as BoFRIa-5 and is also present in the European winter 
type and Chinese semi-winter type parental lines of the 
TN mapping population [22]. 

Functional analysis of BoFRIa alleles in A thaliana 

To ascertain if the two most common BoFRIa alleles 
conferred any functional differences we undertook 
transformation experiments. The coding and 3'UTR 
sequences from the BoFRIa-1 and BoFRIa-4 alleles were 
used to replace the AtFRI coding and 3' UTR sequences 
in an A. thaliana genomic clone. By retaining common 
regulatory sequences in the 5' region from the AtFRI 
gene we hoped to normalise expression and thus focus 
on the structural differences between the two Brassica 
proteins. These constructs were transformed into the 
rapid-cycling A. thaliana accession Columbia (Col-0). 
Col-0 carries a loss-of-function mutation within AtFRI, 
but has a functional FLC so these experiments would 
determine if BoFRIa could complement the fri mutation 
in Col-0, and induce late flowering. Both BoFRIa alleles 
complemented the loss-of-function mutation with > 100 
primary (Tl) transformed plants containing each of the 
BoFRIa alleles flowering very late compared to Col-0 
plants and surprisingly also later than Col-0 transformed 
with a functional AtFRI (Figure 5). 

To investigate the functionality of BoFRIa-1 and 
BoFRIa-4 alleles under different environmental conditions 
five transformants carrying each allele were analysed in 
the next (T2) generation. Flowering time was analysed as 
days-to-flower and total leaf number in plants that had no 
vernalization, or had experienced two or four weeks verna- 
lization, at either 5°C or 10°C (Figure 6). In all treatments, 
except two weeks at 10°C (2W10°C), plants with either 
BoFRIa allele flowered as late as those carrying AtFRI. At 
2W10°C plants carrying either of the BoFRIa alleles flow- 
ered later than AtFRI. Figure 6 also shows that plants 
undergoing a vernalization treatment at 10°C compared to 
5°C continue to grow and initiate leaves at a faster rate. 



Thus, when considering total leaf number as a measure of 
flowering time it appears that only 4W5°C was an effective 
vernalization treatment (Figure 6). 

Expression of the coding sequences of the two BoFRIa 
alleles under the AtFRI 5' regulatory sequences showed 
that both alleles can produce functionally equivalent pro- 
teins that may, under some environmental conditions, be 
even stronger with respect to flowering time effects than 
the endogenous A. thaliana protein (Figure 6). In con- 
trast, two A. lyrata FRI alleles conferred a quantitative 
difference in flowering time by both association and 
transgenic studies [16]. The maintenance of both A. lyr- 
ata alleles at intermediate frequencies in natural popula- 
tions suggests they are differentially selected in different 
environments. If the BoFRIa alleles do underlie flowering 
time QTL then there must be expression differences 
between the two genes to account for the difference in 
flowering time. Both these genes could be expressed in a 
very different pattern to AtFRI as the rearrangement that 
moved it to chromosome 4 resulted in completely differ- 
ent 5' sequences less than 1 kb upstream of the transcrip- 
tion start site and places it in a very different chromatin 
context since it is now 200-300 kb downstream of the 
heterochromatic nucleolar organizer region NOR4 [61]. 

Conclusions 

Knowledge of B. oleracea FRI gene number, functionality 
and map position now puts us in a strong position to 
undertake an extensive investigation into the contribution 
of allelic variation at FRI to flowering, vernalization and life 
history behaviours. Differences in life history between A. 
lyrata and A. thaliana such as outcrossing versus selfing 
and a perennial compared to annual habit may result in a 
requirement for some level of FRI functionality in A. lyrata 
that is optional in A. thaliana [16]. B. oleracea, like A. lyr- 
ata, is a largely outcrossing species and some wild B. olera- 
cea, thought to be the progenitor of the modern crop 
plants, have been reported to keep flowering for up to 20 
years [62]. Our analysis of BoFRIa suggests that only a 
small number of functional BoFRIa alleles are captured 
within the cultivated B. oleracea germplasm. To date we 
have found no evidence for loss-of-function mutations that 
are frequent in AtFRI. Further analysis of the 5' and 3' reg- 
ulatory regions of BoFRI is now underway. The proximity 
of BoFRI to BoFLC, BoFY and BoCO opens up new ques- 
tions of how this may influence flowering behaviour. It will 
be particularly important to be in a position to select speci- 
fic alleles in breeding programmes to allow us to enhance 
robustness against increasing climate variability. 

Methods 

Cloning BoFRI genes 

The JBo BAC library was hybridised with the AtFRI 
genomic clone, originally from the accession Stockholm 
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Caulifbwer E03 
Caulifbwer E05 
Autumn Caulifbwer D10 
Autumn Cauliflower E02 
Autumn Cauliflower Dll 
Autumn Cauliflower D08 
Autumn Cauliflower D12 
Autumn Caulifbwer E01 
Summer Cauliflower E07 
Sumner Cauliflower E08 
Green Cauliflower F05 
Green Cauliflower F06 
Romanesco Caulifbwer F07 
Winter Cauliflower Ell 
Winter Cauliflower E12 
Winter Caulifbwer F02 
Ftoscoff Winter Cauliflower E06 
Broccoli A07 
Broccoli All 
Broccoli B01 
Broccoli B02 
Purple Head Broccoli B08 
Purple Sprouting Broccoli B09 
Feather Leaf Broccoli B06 
Cabbrese B03 
Cafcabrese B05 
Cabbage G07 
Cabbage Cll 
Cabbage D01 



Savoy Cabbage D05 
White Cabbage D07 
Summer Cabbage D06 
Fodder Cabbage D02 
Brussels Sprout Bll 
Brussels Sprout B12 
Brussels Sprout 001 
Brussels Sprout O02 
Brussels Sprout 003 
Brussels Sprout 004 
Brussels Sprout 005 
Brussels Sprout 006 
KaleG03 
KateG04 
Ornamental Kale G08 
Borecole Kale F10 
Fodder Black Kab F12 
Kohl rabi G10 
Kohl rabi Gil 
Kohl rabi G12 
Chinese Kale A02 
Chinese White Kale A05 
Tronchuda Cabbage H02 
Tronchuda Cabbage H05 
Wild Cabbage H08 



Conserved domain 



MIRRREAOATVETVPTNIETTIEOSNDPOFIKSIVDL^ 
MIRRREAOATVETVPTNIETTIEOSNDPOFIKSIVDL^ 

MIRRREAOATVETVPTNTETTTEOSNDPOFT KSTVDT TAT AAAVNAFKRRYDFIOSHMDYTENATDSNT KTNGTVET AAVSPPPKDASGETATAT VCOSPPKEKSEAERT ,CESMCSKF,T RRYMFVNTSERAKT TFFT PGAT KT AKDPAKFVT DCTGKFYT OGRKAFANDSPATTARKVST ,1 ,VT ECYT ,T iTFDPEGEKKOVGSSVKDEAEAAA 

MIRRREAOATVETVPTNIETTIEOSNDPOFTJ<SIVDLTALAAAVNAFKRRYDETJ3SHy^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFl^SIVDLTALAAAVN^ 

MIT^REAOATVETVPTNIETTIEOSNDPOFIj<SIVDLTALAAAVNAFT^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFTJ(SIVDLTALAAAVNAFKRRYDETJDSH^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFTJKSIVDLTO^ 

MIRRREAOATVETVPTOTETTTEOSNDPOFTJ<SIVDLTALAAAVNAFKRRYDETX)SH^ 
MIFRREAOAT\/ETVPTNIETTIEOSNDPOFI^SIVDLTAIAAAVNAFK^ 

MTRRREAOATVETVPTNTETTTEOSNDPOFT KSTVDT TAT ,AAAVNAFKRRYDET DSHMDYTENATDSNT KTNGTWTAAVSPPPKDASGFTATATVOOSPPKEKSEAFRT .CESMCSKELRRYMFVNTSERAKT TFFT PGAT ,KT AKDPAKFVT DCTGKFYT OGRKAFANDSPATTARKVST I VT ECYT T iTFDPEGEKKOVGSSVKDEAEAAA 

MIRRREAOATVETVPTOIETTIEOSNDPOFTKSIVDLTALAAAV^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFI^SIVDLTAIAAAVNAFKRRYDEL^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFLKSIVDLTAIAAAVNAF^^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFIj<SIVDLTAIAAAVNAF^ 

MIRRREAOATVETVPTNIETTIEOSNDPOT^TKSIVDLTALAAAW^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFIJKSIVDLTALAAAVNAT^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFTJ(SIVDLTATJ\AAVNATfK^ 

MTRRREAOATVETVPTNTETTTEOSNDPOFT KSTVDT ,TAT AAAVNAFKRRYDET DSHMDYTENATDSNT KTNGTVET AAVSPPPKDASGETATATVCOSPPKEKSEAERT GESMCSKFT RRYMFVNT SERAKT TFFT PGAT ,KT AKDPAKFVT .DCTGKFYT OGRKAFANDSPATTARKVST ,1 VT ECYT T ,TFDPFGEKKOVGSSVKDEAEAAA 

MIRRREAOATVETVPTOIETTIEOSNDPOFTJ<SIVDLTALAAAVNAFKT^ 

MIFRREA0AT\/ETVPTOIETTIE0SNDP0FIK5IVDLTALA^ 

MTT^REAOATVETVPTNIETTTEOSNDPOFTJKSIVDLTALAAAVNAFT<RRYDETjr)Sr^ 

MIT4RREAOATVETVPTNIETTIEOSNDPOFIj<SIVDLTAIAAAVNAET^ 

MIPRT^AOATVETVPTNIETTIEOSNDPOFT^SIVDLTALAAAV 

MIRRREAOaIvET TIEOSNDPOFIj^SIVDLTAIAAAVNAFKT^YDElJDSH^ TATAIAYOSPPKEKSEAERITESbCSK^^ 

MIT^RREAOATA/ETVPTNIETTIEOSNDPOFI^SIVDLTAIAAAVNAET^ 

MTRRREAOaIvFT TTEOSNDPOFT KSTVDT TAT AAAVNAFKRRYDET DSHMDYTENATDSNT KTNGtIeTAAASPPSPNK TATATAYOSPPKEKSEAERT GESMCSKETRRYMFVNTSERAKT TFFT PGAT ,KT AKDPAKFVT .DCTGKFYT OGRKAFANDSPATTARKVST ,1 VT ECYT T iTFDPEGEKKOVGSSVKDEAEAAA 

MIRRREAOaIvET TIEOSNDPOFT^SIVDLTAIAAAVNAFKT«YDEIX>SHMDYIENAm TATAIAYOSPPKEKSEAERLCESMCSKETJ^YMFVNISERAKLIEELPGALKI^ 

MIRRREAOAIVET TIEOSNDPOFl^SIVDLTAUVAAVNAFKRRYDELOSHMDYIENAIDSNLKTNGliEIAAASPPSPNK TATAIAYOSPPKEKSEAERLCESMCSKELRRYMFVNISERAKLITiELPGALKTJ^ 

MIRRREAOATVETVPTNIETTIEOSNDPOFLKSIVDLTAIJ\AAVNAT^ 
MIRRREAOATVET TIEOSNDPOFIJ<SIVDLTAIAAAVNAFT<T«YDElJ3S 



mtrrreaoa vet- 
mirrreaoa vet- 
mtrrreaoa vet— 
mirrreaoa vet— 
mirrreaoa vet— 
mtrrreaoa vet— 
mirrreaoa vet- 
mirrreaoa vet- 
mtrrreaoa vet— 
mirrreaoa vet— 
mirrreaoaIvet— 



— TIEOSNDroFTJKSTVDLTATJ\AAVNAFT<RRYDETJDSHMDYIENAIDSNLKTNGI EIAAASPPSPNK 

-TIEOSNDPOFTJ<SIVDLTATJVAAVNAFKRRYDETJDSHMDYIENAIDSNLKTO^ EIAAASPPSPNK- 
-TTEOSNDPOFT KSTVDT TAT AAAVNAFKRRYDET DSHMDYTENATDSNT KTNGT 

--TIEOSNDPOFU^IVDLTALAAAVNAFKRR^^ EIAAASPPSPNK- 

-TIEOSNDPOFT^SIVDLTATJ\AAVNAFT<RRYDETJDSHMDYIENAIDSNLKTNGI eiaa? 



P K— 

PNK— 
PNK- 
PNK- 
PNK— 
S P K— 
EIAAASPPSPNK— 




:TVP™iETTIEOSNDPOFTJ<SIVDLTATJVAAVNAFT<T«YDEmSHMDYTENAIDSNLKTOGIV 
■TIEOSNDPOFT^SIVDLTATJVAAVNAFKT«YDETX)SHMDYIENAIDSNLKTNGll 
EOSNDPOFT KS TVDT TAT .AAAVNAFKRRYDET DSHMDYTENATDSNT KTNGtI 



TIEOSNDPOFTJ(SIVDLTAIAAAVNAFKP^YDETJDSHMDYIENAIDSNLKTNGli 



mtrrreaoatvftvptntettteosndpoft kstvdt tat ,aaavnat^rydftoshmdytenatdsntktngtvf," 
mirrreaoaIvet tieosndpoftj^sivdltatj\aavnafkrrydetjdshmdyienaidsnlktngibe: 



TATAIAYOSPPKEKSEAERLmSiynSKETJ^YMFvTSIISERAKLIEELK 

tataiayosppki^seaeri££smcsketj3rymfvniserm 

t aaasppspnk tatatayosppkekseaert gesmcskft rrymfvnt serakt tfft pgat ,kt akdpakfvt .dctgkfyt ogrkafandspattarkvst i ,vt ecyt ,t itfdpegekkovgssvkdeaeaaa 

tataiayosppkekseaerij^esmcskeijrrymfvniseraklit^ 
tataiayosppjkt^seaerixesmcskeh^ymfvnisfj^ak^ 

tatatayosppkekseaert gesmcskft rrymfvnt serakt tfft pgat kt akdpakfvt .dctgkfyt ogrkafandspattarkvst .1 vt ecyt i itfdpegekkovgssvkdeaeaaa 
tataiayosppki^seaerlcesmcskeijrrymfvnisepakl^ 
tataiayosppkekseaerlcesmcskei^ymt^iserakxieelpgalklakdpakfvltx 
tataiayosppkekseaertjtttismoskeijrrymfvniseraklieel^ 
tataiayosppkekseaerlffisncsketj^yme^iserak^ 
tataiayosppkekseaerlcesmcskeijrrymfvniseraklieelpga 
iaavspppkdasgetataiacosppkokseaerlcesmcskeijrrymfvniserat^ 
iaaasppspnk tataiayospplktitkseaertxesmcskeij^ymfvniser 

t aaasppspnk tatatayosppkekseaert gesmcskft rrymfvntserakt tfft pgat ,kt akdpakfvt .dctgkfyt ogrkafandspattarkvst i vt ecyt t itfdpegekkovgssvkdeaeaaa 

iaaasppspnk tataiayosppkekseaet^ijcesmcskeijrrymfvniserak^ 

iaaasppspnk tataiayosppikekseaerlcesmcskeiroymfvniserai^ 

iaaasppspnk — tataiayosppkekseaerigesmcsketj^ymfvnisefaklieel^^ 

iaaasppspnk tataiacospfkekseaerlffisncskei^ 

iaavspppkdasgetataiacosppkokseaerlcesmcskeit^rymfvniset^ 

iaaasppspnk tataiayosppkekseaerioesmcskeijrrymfvniseraklit^l^ 

iaaasppspnk tataiayosppkt^seaerixesmcskeij^ymfvnisfj^ak^ 

t aavspppkdasgftatat acosppkokseaert gesmcskft rrymfvntserakt tfft pgat kt akdpakfvt dctgkfyt ogrkafandspattarkvst ,t vt ecyt t itfdpegekkovgssvkdeaeaaa 
iaaasppspnk tataiayosppkekseaerlcesmcskeij^ymevniseraklieelpgai^^ 



Figure 4 Comparison of first 160 amino acids of BoFRIa from 55 BolDFS genotypes. Protein comparison of BoFRIa from cultivated genotypes of B. oleracea listed by crop type. The plate 
co-ordinates refer to those listed in Table 1. The conserved region of 37 amino acids that defines FRI from other members of the FRI superfamily is delineated by the horizontal black bar. 
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Figure 5 Distribution of flowering times in T1 transformants 
carrying BoFRIa-1 and BoFRIa-4 alleles. Histogram of the 
flowering time of T1 lines transformed with BoFRIa-1 and BoFRIo-4 
measured as days to flower. The flowering time of wild type Col-0 
plants and Col-0 transformed with AtFRI are indicated. 



[8] and seven BACs identified, six having identical 
restriction patterns and one different. Purified genomic 
DNA was prepared (Qiagen Maxi Prep Kit) from two of 
these BACs (72123 and 88G16) and used to generate 
shotgun libraries (TOPOr Shotgun Kit) of 1-2 kb frag- 
ments, in the pCR R 4Blunt-TOPO R vector, giving 6-fold 
coverage. Colonies from these libraries were gridded 
onto nylon membrane (HyBond-N+) and hybridised to 
three probes generated from AtFRI (the 5' region, exon 
1, and the 3' end of exon 3 and 3' UTR). BAC sub- 
clones were identified with each of these probes. 
Sequence analysis confirmed that the two BACs carried 
different Brassica paralogues. 

Mapping BoFRI loci in B. oleracea mapping populations 
Genetic mapping 

Mapping filters of the A12DHdxGDDH33 mapping popu- 
lation were produced and hybridized with a conserved 
BoFRI probe as described in [36,63]. A 900 bp conserved 
region from exons 2 and 3 from BoFRIb was amplified 
from A12DHd genomic DNA with primers J2NG_F1 (5' 
AAGTATCAAGCGTGGAAAGCA 3') and J2NG_R1 (5' 
GTTACGAGGAGACCTGTGATT 3') and used to probe 
both the A12DHdxGDDH33 and NxG mapping filters 
(supplied by Graham Teakle, WHRI). Linkage analysis to 
map the BoFRIa locus was performed using Joinmap 3.0 
[64] with the mapping data provided at BrassicaDB http:// 
brassica.bbsrc.ac.uk/BrassicaDB/. 
Fluorescence in situ hybridisation (FISH) 
FISH was performed on chromosome spreads from the 
A12DHd genotype of B. oleracea using methods described 
in [65]. The chromosomes are now named according to 
their corresponding linkage group. JBo72I23 was applied 
to meiotic pachytene spreads together with JBo62M08, a 
BAC which is associated with the RFLP marker pN22 on 
C3 at 42 cM and previously assigned to chromosome C3 
by FISH. JB088GI6 {BoFRIb) was applied to mitotic meta- 
phase spreads together with BAC BoB061G14, which 
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Figure 6 Functional analysis of the two most common BoFRIa 
alleles. Average flowering time of T2 families transformed with the 
two major BoFRIo alleles BoFRIa-1 and BoFRIo-4 and compared to 
Col-0 transformed with the AtFRI allele and Col-0 wild type, (a) 
Flowering time measured as days-to-flower. The error bars show 
95% confidence intervals, (b) Flowering time measured as final leaf 
number. The error bars show 95% confidence intervals. Segregating 
progeny with and without the transgene are indicated by + and - 
respectively. 



hybridizes to pericentromeric heterochromatin of six pairs 
of chromosomes, and a 45S rDNA probe from clone 
pTa71 [66], EMBLX07841. The chromosome pair to 
which JB088GI6 hybridized lacked signals from the other 
probes and had morphology suggestive of C9. Therefore, 
slides were reprobed with JBo32J18, a BAC associated 
with BoFLCl which has been mapped to a region between 
pN47E4NM (87 cM) and pN3El (103 cM) on C9 [31,67] 
and confirmed to be on C9 by FISH (unpublished). Two 
further BACs showing the same restriction pattern as 
JB088GI6 were applied separately with JB088GI6 to 
pachytene spreads. 

Sequencing BoFRIa in BolDFS 

The B. oleracea diversity foundation set (BolDFS) is a core 
collection of lines that represent the genetic variation 
across the morphologically diverse crops of this species 
http://www.brassica.info/resource/plants/diversity_sets. 
php. DNA was isolated using the DNeasy 96 Plant Kit 
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(Qiagen) and amplified using the GenomiPhi whole gen- 
ome amplification kit (GE healthcare). A 650 bp fragment 
of BoFRIa was amplified from genomiphied DNA of 55 
genotypes of the BolDFS by PCR with primers YWFRI_F 
(5'CGCACATCGTCCATCAACAAG 3') and FRIJ1_R2 
(5'ATCCTTCACCCACCAGCCT 3') using AMPLITAQ 
GOLD TAQ DNA Polymerase (Life Technologies Ltd 
(Invitrogen Division)). Sequence analysis was conducted 
using AlignX in Vector NTI (Invitrogen). 

Functional analysis of BoFRIa alleles 

Plasmid pFRIg (in pBluescript-KS+, Stratagene) was 
mutagenised to introduce a BamHl site immediately 5' 
of the ATG (plasmid pFRIg-B). Digestion of pFRIg-B 
with BamHl plus Clal allowed removal of the AtFRI 
coding sequences, leaving the 5' region of AtFRI. A 4.3 
kb fragment containing BoFRIa was isolated from geno- 
mic DNA of lines El and E8 by PCR with primers 
BoFRIl_ Bam_ATG (S'CTTCCGCGGATCC- 
CATGGCCGTCCGTAACS') and BoFRIl_R2_ClaI 
(5'CAGAGATCGATC TCGAGAAAGGTAGCTGTTT 
3'), using Pfullltra II Fusion HS DNA Polymerase (Agi- 
lent Technologies) and sequenced. PCR products were 
digested with BamHl plus Clal and the purified frag- 
ments ligated into BamHl plus CM-digested pFRIg-B to 
give final constructs containing the 5'UTR of AtFRI with 
the coding and 3' UTR sequences of BoFRIa-1 (in 
BoFRIa-1) or BoFRIa-4 (in BoFRIa-4). pFRIg-B was used 
as the A. thaliana FRI control. The final constructs 
were ligated into binary vector pSLJ755I6 (a gift from 
Prof. Jonathan Jones, http://www.tsl.ac.uk/research/ 
jonathan-jones/plasmids.htm), on an EcoRl plus Xhol 
fragment (pFRIg-B) or EcoRl plus Clal fragments (from 
BoFRIa-1 and BoFRIa-4). The constructs were trans- 
ferred into Agrobacterium by triparental mating [68] 
and transformed into A. thaliana accession Col-0 by a 
floral dipping method (modified from [69]). Tl transfer- 
mants were isolated by selection for Basta™ resistance. 
T2 seed were collected and flowering time determined 
by days-to-flower excluding the period of vernalization 
treatment and final leaf number at flowering. 

Plant growth 

T2 A. thaliana seeds were sown on 'Arabidopsis mix' 
(Scotts® Levington F2 8.75 1 bags, 100 1 of grit, 200 g of 
Imidasect® 5 gr.) in plastic pots (7 cm x 7 cm) and strati- 
fied in a vernalization chamber at 5°C with an 8 h photo- 
period and constant humidity for 3 days. Pots were moved 
to a naturally lit long day glasshouse for 7 days in May 
2010 to allow germination and pre -growth. Seedlings not 
receiving a vernalization treatment remained in the glass- 
house; seedlings to be vernalization treated were trans- 
ferred back to a vernalization room or controlled 
environment cabinet (Snijder Economic Deluxe) for a 



treatment of either two or four weeks at 5°C or 10°C. 
After vernalization, 20 plants per line were transplanted 
into trays with 40 cells of 2 cm x 2 cm and returned to 
the glasshouse. Trays were moved regularly to random 
positions to prevent any positional effects on plant growth. 
Flowering time was recorded as either total leaf number 
(rosette leaves plus cauline leaves at flowering) or bolting 
time; bolting time was scored as the number of days-to- 
flowering determined when the inflorescence stem was 
3 cm tall. 

Abbreviations 

QTL: Quantitative trait loci; BAC: Bacterial artificial chromosome; ORF: open 
reading frame; FISH: fluorescence in situ hybridization. 
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